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ABSTRACT 

We use a high-resolution dissipationless simulation to study the evolution of the dark 
matter and halo distributions in a spatially flat cosmological model dominated by a 
cosmological constant A and cold dark matter (ACDM). In order to quantify the evo- 
lution of structure, we calculate the Minkowski functionals of the halos and the dark 
matter component at various redshifts. A comparison of Minkowski functionals and 
the more standard correlation function analysis shows that the Minkowski functionals 
contain information about correlation functions of arbitrary order. While little evolu- 
tion of the Minkowski functionals of halos between z = A and 2; = is observed, we 
find that the Minkowski functionals of the dark matter evolve rapidly with time. The 
difference of the Minkowski functionals of halos and dark matter can be interpreted 
as a scale dependent bias. This implies that scale-dependent bias is a property of not 
only the two-point halo correlation function, but also of correlation functions of higher 
order. 

Key vifords: Methods: statistical; Galaxies: general; Cosmology: theory; dark matter; 
large-scale structure of Universe. 
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1 INTRODUCTION 

It is generally believed that dark matter (DM) constitutes a 
large fraction of the mass in the Universe. Therefore, it sig- 
nificantly affects both the process of galaxy formation, and 
the large-scale distribution of galaxies. The gravitational 
domination of DM on the scale of the galactic virial radius 
implies that collisionless simulations can be used to study 
the formation of the DM halos of galaxies. 

The standard flat cosmological scenarios with cold dark mat- 
ter alone cannot explain the structure formation both on 
small and very large scales. Recently, scenarios with a non- 
zero vacuum energy, quantified by the cosmological constant 
A, have been proved to be very successful in describing most 
of the observational data at both low and high redshifts. 
Here we study the evolution of the dark matter and halo 



distribution, a variety of statistical methods have been devel- 
oped. The hierarchy of correlation functions, counts-in-cells 
or the power spectrum are examples of low-order statistics. 
While this may be sufficient to characterize the large-scale 
structure of the Universe, gravitational clustering quickly 
leads into the nonlinear regime where strongly non-Gaussian 
structures evolve. More information can then be expected 
from measures that provide a handle on the global mor- 
phology of the large-scale structure. One example of such 
measures are the Minkowski functionals. 

Minkow ski functionals have b een introduced into cosmol- 
ogy by 



Mecke et al. (1994), They provide a geometri- 



distribi tions in a spatially fiat cosmological model domi- 
nated by the cosmological constant and cold dark matter, 

acdmI 



cal and topological descripti on of a point distributio n (for 
a detailed introduction, see Schmalzing et al. 199q ). Re- 
cently Minkows ki functionals have been applied to clus- 
ters of galaxies (iKerscher et al. 1997), IRAS galaxies (Ker 



scher et al. 1998D , an d anisotropy maps of the cosmic mi- 
crowave background ([Winitzki fc Kosowsky 1997; Schmalz- 



ing fc Gorski 199^ ; Novikov et al. 1998) 



In order to assess the properties of the dark matter or halo 
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Our article is organized as follows. In section Q we summarize 
the properties of our ACDM simulation, while in section p^ we 
briefly introduce the method of Minkowski functional analy- 
sis. In the following section W we apply Minkowski function- 
als to describe the distribution of dark matter particles and 
halos. Finally, in section H we summarize obtained results 
and draw our conclusions. 



2 THE SIMULATION 

We investigate a spatially flat model with a cosmological 
constant (ACDM with Q.o = 0.3, Q.a = 0.7, as = 1.0, and 
h = 0.7). These parameters allow for a power spectrum nor- 
malization in accord with both th e four year COBE DMR 
observations ( Bunn fc White 1997) and the ob served abun- 
dance of galajcy clusters (Viana & Liddle 1996). The age of 
the universe in this model is ~ 13.5 Gyrs. 



Colin et al. (1999). The BDM code is available from authors 



In order to study the statistical properties of halos in a cos- 
mological environment, the simulation box should be sufB- 
ciently large; we use a cube of 60/i~^Mpc side length. On 
the other hand, to assure halo survival in clusters, the force 
resolution should be < 1 — 3 fc~^kpc and the mass resolu- 
tion shou ld be < lO^/i'^M© ( [Moore et al. 1996[ [Klypin et 
al. 1999 ). Our simulations we re done using the A daptive 
Refinement Tree (ART) code (Kravtsov et al. 1997) with a 



dynamical range of 32,000 in high density regions. The simu- 
lation followed evolution of 256'^ cold dark matter particles, 
which give particle mass of 1.1 x 10 h~^ Mq. The mass res- 
olution is therefore sufficient to identify galaxy-size halos of 
mass Af > 5 - 10 X lO^^^-^M©. 

Identification of halos in dense environments and reconstruc- 
tion of their evolution is a challenge. A variety of situa- 
tions that are frequently found in the real Universe are dif- 
ficult to identify in simulations. Typically, problems arise if 
a small satellite is bound to a larger galaxy; famous exam- 
ples are the M51 system, or the LMC and the Milky Way. 
Cases when many small gravitationally self-bound halos are 
moving within a large region of virial overdensity, such as 
clusters and groups of galaxies, can also cause trouble. We 
have developed two algorithms to deal with such situations, 
namely the hierarchical friends-of-friends (HFOF) pr ocedure 
and the bo und density maxima (BDM) algorithm (Klypin 
et al. 1999). Both yield consistent results, i.e. they basically 



identify the same halos in a DM simulation. 

In this study, we use the BDM algorithm to construct halo 
catalogs at different redshifts. The main idea of the algo- 
rithm is to find positions of local maxima in the density 
field smoothed on a certain scale and to apply physically 
motivated criteria to test whether the identified site corre- 
sponds to a gravitationally bound halo. The density max- 
ima are found iteratively for a series of smoothing scales, the 
smallest scale being close to the peak resolution of the simu- 
lation. To test whether a density maximum corresponds to a 
real halo, we construct density and velocity dispersion pro- 
files centered on the maximum and iteratively remove un- 
bound particles. This procedure eliminates unbound halos 
from the catalog, while also cleaning bound halos off some 
high-velocity background particle s. The detailed desc ription 
of the BDM algorithm is given in (Klypin et al. 1999) and in 



upon request. 

Halos are usually characterized by their virial mass. Unfor- 
tunately, it is difficult to obtain a meaningful definition of 
this quantity for halos inside larger halos, so we avoid the 
problem by determining the maximum circular velocity Ucirc 
of each halo. The circular velocity is defined as the maxi- 
mum of ^GM{< R)/R, where M(< R) denotes the mass 
contained inside a sphere of radius R around the halo cen- 
ter, and G is the gravitational constant. Vdic is more useful 
both observationally and numerically. Moreover, it can be 
measured more easily and more accurately than mass. For 
isolated halos the maximum circular velocity and the virial 
radius are, of course, directly related and are therefore equiv- 
alent. 

The output of the halo finding algorithm is primarily de- 
termined by the assumed lower mass threshold and also de- 
pends weakly on the assumed maximum halo radius. With 
the threshold of 10^°/i"^A'/q (10 particles) at 2 = the algo- 
rithm identifies about 17,000 halos in our simulation. In or- 
der to conduct a study of halo statistics one needs a complete 
halo sample that is not affected by the numerical details of 
the halo finding procedure. We have tested the completeness 
of the halo samples w ith the differential velocity functions 
([Gottlober et al. 1998). Depending on redshift one can de- 



fine a threshold of the circular velocities so that for circular 
velocities larger than this threshold the halo samples do not 
depend on the numerical parameters of the halo finder. It 
has been verified that at 2 = the sample is complete for 
all halos with Ucirc ^ 85km/s, while for 2 = 1 completeness 
is achieved only with iicirc > lOOkm/s. 



3 MINKOWSKI FUNCTIONALS 

Let us consider an object in three-dimensional Euclidean 
space. Its morphology can be characterized by the four 
Minkowski functionals, namely the volume V, the surface 
area A, the integral mean cur vature H , and the Euler char- 
acteristic X- Hadwiger (1957) has shown that the Minkowski 
functionals supply a complete and unique description of the 
global morphology. As geometrical characteristics of struc- 
ture, the Minkowski functionals combine both the advantage 
of a intuitive geometrical interpretation and the advantage 
of delivering a quantitative measure. 

In the following, we look upon the distribution of the dark 
matter particles and the halos as realizations of point pro- 
cesses. To investigate the distribution of these points we dec- 
orate each of them with a ball of radius r. Starting from 
the trivial situation at r = 0, where all points are isolated 
from each other, we increase r, thereby creating connections 
between neighboring balls. This procedure results in a com- 
plex pattern of intersecting balls the properties of which can 
be studied using Minkowski functionals. To be precise, we 
calculate for a set of points {xj} decorated by balls B of ra- 
dius r the Minkowski functionals Vp ( |J Bx ) of the union 
set of all balls. Varying the radius introduces a diagnostic 
parameter into our investigation. Since the set of points is 
generated by a reasonably well-behaved point process (see 
Weil 1983 for a rigorous discussion of requirements) , we can 
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define an ensemble average of the Minkowski functionals per 
unit volume. We call these quantities n^. 

Like other morphological estimators, the Minkowski func- 
tionals of a point sample depend on its number density. In 
the case of a Poisson process one may derive scaling rela- 
tions with the number density and ra dius using the analyti- 
cal results of (Mecke & Wagner 1991) and the homogeneity 
property of the Minkowski functionals. For general distri- 
butions, however, scaling relations are not readily available 
since the scaling properties depend, even without correla- 
tions, on the dimensionality of the support of the point pro- 
cess. This is why we always take care to select the same 
number of points from our samples before we calculate the 
Minkowski functionals. However the total number of halos 
above a certain minimum circular velocity changes with red- 
shift due to merging. Therefore, we proceed as follows: We 
assume a threshold of lllkm/s for the circular velocity Vdic, 
and find 5575, 9644, 5869 halos at redshifts 4, 1, and 0, re- 
spectively. From these halos we randomly select at each mo- 
ment 5000 halos, and compute their Minkowski functionals. 
We do this several times with different random sets to de- 
rive a scatter. To compare the distribution of halos with that 
of dark matter we also select 5000 dark matter particles at 
random. 

All the point samples we study live in cubic boxes with pe- 
riodic boundary conditions. Her ice correction of edge effec ts 



as performed, for example, by gchmalzing et al. (1996) in 
the analysis of real catalogues, is not an issue here. The 
numerical computation of the Minkowski functionals is per- 
formed usi ng an impleme ntation of the partition formula 
derived by Federer (1959). The formula provides a general- 
ized concept of surface integration, suitable for taking into 
account the corners and edges encounte red in our union se t 
of balls (for a detailed description, see 



Mecke et al. 1994) 



The present implementation offers the advantages of scaling 
linearly with the number of points, and requiring only mod- 
erate computational resources; a sample of 5000 points can 
be fully analyzed in about one hour on an average worksta- 
tion. Together with a parallelized version, the source of the 
code used in this work is available from the authors upon 
request. 



4 ANALYSIS 

4.1 Evolution of two-point correlation function 
and Minkowski functionals 

Recently, the evolution of clustering of dark matter halos 
has been investigated for several cosm ological models using 



numerica l dissipation l ess simulations ( Brainerd fc Villum- 



sen 1994 



Bagla 1998| ; [Colin et al. 1999| ; [Kravtsov fc Klypin 



199£; Ma 199£), simulations that include both dissipation 



less dark matter and dissipative baryonic component s 



et al. 1999; Blanton et al. 



1998, 



Blanton et al. 




Ostriker 199£ ; Jenkins et al. 1998D, and "hybrid" studies in 



which dissipationless simulations are compl emented with a 
semi-analytical model of galaxy formation ([Roukema et al 



1997|; iKauffmann et al. 1999a|; [Kaufi'mann et al. 1999bt |Di| 



aferio et al. 1999| ; [Baugh et al. 1999| ; [Benson et al. 1999| ; [Ko, 
latt et al. 19991). AH of these studies, though very diverse in 
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Figure 1. Correlation functions of the iialo samples (fcirc > 
lllkm/s) and the dark matter distributions at z = 0, 1, 4. 

their methods, qualitatively agree on one important result: 
bias is expected to be non-linear, to depend on the proper- 
ties of the DM halos and the galaxies they host, and to be 
a (generally non-monotonic) function of cosmic time. 

In particular, Colin et al. (1999) discuss the evolution of the 



real space two-point correlation function in different cos- 
mological models. They find that in all models and at all 
epochs, the correlation function of halos can be reasonably 
approximated by a power-law. In Figure [ll we compare the 
correlation function of the dark matter particles in real space 
with the correlation function of halos with circular veloci- 
ties «circ > lllkm/s. In agreement with previous studies, 
we see almost no evolution of the halo correlation function 
between redshifts 4 and 0, whereas by comparison the cor- 
relation function of the dark matter particles increases with 
time. 

The four Minkowski functionals for halos and dark matter 
particles at various epochs are shown in Figures H and pi 
respectively. The scatter between different subsamples is 
marginally larger than the line width, so for clarity we only 
plot the mean values. One should keep in mind that due 
to the density dependence of the Minkowski functionals the 
diagnostic parameter "radius" cannot be directly associated 
with a scale. Instead, the general features described below 
occur within the radius range between zero and a few times 
the average distance for any point distribution studied. It is 
the location and relative strength of these features that ac- 
counts for the discriminative power of the Minkowski func- 
tionals. Quantitative measurements may be possible using 
reference distributions; the simplest such case is the Poisson 
process displayed in all Figures. 

The zeroth Minkowski functional is nothing but one less 
the void probability function (VPF), which is defined as 
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Figure 2. Minkowski functionals of the halo samples (fcirc > lllkm/s) samples at 2 = (solid), z = 1 (dashed) and 2 = 4 (dotted). 
Note that the results for 2 = 1 and 2 = 4 coincide almost everywhere, so the lines appear dash-dotted. As a standard of reference, we 
show the analytically known values for a Poisson process of the same intensity (thin solid line). 



the probability of finding no objects within rando mly placed 
balls (for a recent study of the VPF of halos, see Ghigna et 
al. 199n) . Voids are successively filled as the radius of balls 
increases, so the vo increases monotonically approaching a 
value of one at complete filling of the whole simulation box. 
The first and second Minkowski functionals are related to 



the surface content and the integrated mean curvature of 
the body, and have rarely been applied in cosmology. Both 
quantities start from a value of zero at zero radius, and in- 
crease to a maximum value. While the area only decreases 
from that point, as the union set becomes more connected, 
the integrated mean curvature shows a turnover from pos- 
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Figure 3. Minkowski functionals of the dark matter distribution at z = (solid), z = 1 (dashed) and z = 4 (dotted). As in Figure 
we include a Poisson process (thin solid line). 



itive to negative values, indicating a change from convex 
to concave structures. The third Minkowski functional, the 
Euler characteristic us, is a purely topological quantity and 
counts the number of components less the number of tunnels 
in a structure. At small radii, all balls are isolated from each 
other and v^ is simply equal to the number density of points. 
As connections are established between closely neighboring 



points, the Euler characteristic decreases and eventually falls 
below zero, when an intricate network of tunnels has formed 
in the union set. A second maximum is reached at an already 
fairly large filling factor, when the union set mainly contains 
isolated cavities, which count as components and each con- 
tribute +1 to the Euler characteristic. Eventually the whole 
box is filled and Da becomes zero. Note that the Euler char- 
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acteristic «3 is connected to the genus g via g = 1 — 2«3 . The 
genus has been already widely discussed in the cosmologi- 
cal context using density fields const ructed by fixed-widt h 
smoothing the particle distribution (Gott III et al. 1986). 
The result then depends on both the smoothing scale and 
the density threshold which has been used for the construc- 
tion of the body. 

In Figure ^ we show the Minkowski functionals of 5000 ran- 
domly selected halos at redshifts z = 4, 1, and 0. We do not 
show the scatter resulting from different selection, since the 
standard deviation is of the order of the line width. Only 
little evolution of the clustering of halos can be seen. In 
particular, the Minkowski functionals ai z — 4 and z = 1 
coincide, while the number of halos with Vdvc > lllkm/s 
increases in this time interval by almost a factor of two. If 
one calculates the Minkowski functionals of all halos above 
that threshold of the circular velocity one would see weaker 
clustering at z = 1 than at z — 4. This demonstrates the 
importance of selecting the same number of objects for the 
comparison of Minkowski functionals. 

Already at z = 4 the Minkowski functionals of halos are 
markedly different from a Poisson process. We interpret the 
evolution of the Minkowski functionals in the long time in- 
terval between z = 1 and z = as nonlin ear clustering 
on small scales (see Kravtsov & Klypin 1999 for a detailed 



discussion). The halos tend to move towards existing clus- 
ters of halos, hence the void probability function increases 
on all scales, and both surface and mean integral curvature 
decrease on small scales. 

It is remarkable that the Euler characteristic of the halo 
distribution hardly changes with time. The absence of the 
positive maximum of vz at any time and any radius in the 
halo distribution can be interpreted as the signature of a 
network of filaments or sheets, which leave only few cavities 
in the almost filled union set. Furthermore, the sharp drop of 
the Euler characteristic at small radii points towards strong 
clustering on small scales. These features reveal the pres- 
ence of a network of filamentary structures spanning the 
whole simulation box. We are confident that our result is 
only weakly infiuenced by the comparatively small box size 
of 60/i~^Mpc. After all, two points in the set can only inffu- 
ence each other's contribution to the Minkowski functionals 
if they are less than twice the radius apart. However, com- 
plete filling occurs below 10/i~^Mpc, i.e. almost an order of 
magnitude below the box size. 

In contrast to the modest evolution seen in the halo distribu- 
tion, the clustering of dark matter shown in Figure ^ evolves 
considerably with time. The present day clustering of DM 
appears to be as strong as that of the 2 = halo sample. At 
2 = 4, however, the DM clustering properties still resemble 
those of the Poisson process. The modest differences from 
a purely random distribution refiect that most of the parti- 
cles are still contained in the field, and only a minority has 
formed structures. Even if the sampling density is increased, 
this effect persists. 

Finally, let us compare the Minkowski functionals of the DM 
and halo distributions. The DM distribution evolves continu- 
ously towards higher clustering as one has seen already from 
the analysis of the evolution of the correlation function and 
the power spectrum. At early epochs the Minkowski func- 
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Figure 4. Fo r ha los at z = we compare the asymptotic se- 
ries from Eq. (A.15) (solid line) to its truncation after the first 



non-trivial term, which contains a convolution of the two point 
correlation function. The dashed line is obtained using the halo 
correlation function of Fig. (hh. 
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4r/3 
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(l-df(l-d/2) 
1-d 


2 
3 


1 - d + (Vl - d2/2) arcsind 
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Table 1. Minkowski functionals of a single ball B of radius r 
and of an intersection B n -Bx of two balls separated by a distance 
X = 2rd {0 < d < 1). 



tionals indicate a much stronger clustering of the halos than 
of DM, but there is only a little difference in the morpho- 
logical properties of halo and DM distribution at z = (cf. 

Sect. ka). 



4.2 Comparison to Correlation Function Analysis 

In theory, Minkowski functionals contain information about 
correlation functions of arbitrary order. However, it is not 
clear a priori that this actually results in enhanced informa- 
tion compared to the standard correlation function analy- 
sis. Although the analysis of the previous section revealed a 
wealth of information on the evolution of clustering in our 
simulation, it seems that qualitatively the same informa- 
tion could have been deduced from the correlation function 
alone. In Appendix iM we establish an analytical connec- 
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tion between the hierarchy of correlation functions and the 
Minkowski functionals, and use this relation to conduct a 
direct comparison. 

The Minkowski functionals per unit volume, v^, of a Poisson 
process with intensity p are given by 

pVo(B) 



Vo 

vi = e 
V2 = e 

V3 = e 



1-e 

pVo(B) 



pViiB) 



-pVo(B) 



8 
9 



97r 



pVsiB) - -p^Vr{B)V2{B) + —p'VriB) 



(1) 

where V^(-B) is the /ith Minkowski functional of a ball B 
of radius r. As shown in Appendix W this formula can be 
generalized for arbitrary point processes, if the V^(-B) are 
replaced by a series V^ containing the complet e hie rarchy 
of correlation functions. Truncating Equation (A12) after 
the first two terms we have 



V^ = V^{B) - I / d\^{x)V^{B n B^) ± 



(2) 



In Table Py we summarize the Minkowski functionals Vjj(i3) 
of a single ball B of radius r and of an intersection B n Bx of 
two balls separated by a distance x — 2rd (0 < d < 1). Note, 
that all Minkowski functionals of the empty intersection x > 
2r become zero. 

In Figure W we present the asymptotic series V^ of the 
halo distribution at 2: = (solid line). We co mpar e the full 
asymptotic series obtained by inverting Eq.(All) with its 
truncation after the first non-trivial term (dashed line) ac- 
cording to Eq. (Pp. Here we use the halo two-point corre- 
lation function at z = determined from the simulation as 
presented in Fig. 111. Obviously, the leading terms do not even 
vaguely approximate the full series from Equation (0). 

Integrating Eq. (H) with a standard power-law correlation 
function ^(x) = x ' we obtain the scaling behavior 



l-V^,/V^,{B)<xr 



3-7 



(3) 



for the truncated series. For clarity, we have not included 
these curves in Fig. W but with 7 = 1.7, a reasonable fit to 
the curve calculated from the numerically determined cor- 
relation function can be obtained. This confirms that the 
Minkowski functionals reveal information beyond the corre- 
lation function analysis, and sharply discriminate between 
the halo distribution of our model and a purely Gaussian 
process with the same second-order characteristics. Since 
they reveal clustering properties of a totally different qual- 
ity, Minkowski functionals are an ideal complement to more 
standard methods of large-scale structure analysis. 
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Figure 5. In these plots we compare the ratios 1 — V i_i/Vi_i{B) 
for halos (solid lines) and dark matter particles (dashed lines) at 
redshift z = 0. The areas indicate the standard deviation between 
different subsamples of the point distribution. 
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4.3 Bias 

The comparison of Minkowski functionals on its own pro- 
vides a qualitative understanding of the clustering proper- 
ties of matter in the model. In this section we compare di- 
rectly the clustering properties of DM and halos at z — 0. 
It is well known that the distribution of halos is different 
from the DM distribution. The relationship between these 



© 1999 RAS, MNRAS 666, 0-0 



Figure 6. As in Figure a, these plots show the ratios 1 — 
Vii,/Vi_i{B) for halos (solid lines) and dark matter particles 
(dashed lines), but at redshift z = A. Again, the areas display 
the scatter between subsamples of 5000 points taken from the 
whole particle distribution. 
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two distributions is called bias. Considering the correlation 
functions of halos and DM the ratio between both can be 
interpreted as scale-dependent bias (Fig. m). 

This scale-dependent bias is a generic prediction of hier - 
archical models of structure formation ( polin et al. 1999 ). 
The bias is also time-dependent: it decreases from an ini- 
tially high value at 2; ~ 3 — 7 to a considerably smaller value 
(which can be even less than one) at present. Qualitatively, 
the same behavior has been found in the pow er spectrum 
both in real space (Kravtsov fc Klypin 1999 ) and in red- 
shift space (Gottlober et al. 199J). Since power spectrum 
and correlation function are the Fourier pairs, a bias in one 
of them leads to a bias in the other. Obviously, any scale 
dependence will be different in coordinate and k space. 

In Figures ra and o we compare the functions 1 — l/fi/V^(i3) 
for the dark matter and halos. Our analysis indicates that 
these functions (the solid lines in Figure kl) are much more 
suitable for comparison than the Minkowski functionals 
themselves. For example, 1 — F,j/V'^(-B) = is expected 
in case of a Poisson process. The figure shows these quanti- 
ties for halos (solid lines) and dark matter particles (dashed 
lines). At z = 0, bias is a lot weaker than at z = 4. For 
small radii, an anti-bias is observed, meaning that halos 
appear less clustered than dark matter particles. This does 
not occur at higher redshifts, where halos cluster stronger 
than dark matter particles at all radii. Although the radii 
for which the halo and DM curves coincide are different for 
the four functionals and different from the radius where the 
correlation functions of halos and DM coincide, the qual- 
itative appearance of all Minkowski functionals is similar; 
this reflects the fact that all 1/^ depend on the hierarchy of 
correl ation functions in a very similar way through Equa- 
tion (|A12|). 



of groups and clusters. Indeed, the dynamical friction results 
in mergers of satellite halos with the central cluster halo and 
thereby reduces the ratio of halo overdensity to the ever- 
growing overdensity of dark matter. We refer the reader to 



Kravtsov & Klypin (1999) for an extended discussion and 



The two solid and dashed lines in each of the panels indicate 
the scatter which we have estimated from the statistics for 
different random selections of 5000 objects. At large radii, 
uncertainties increase because the Boolean grain model al- 
ready fills a considerable portion of the whole space, and 
small fluctuations in the pattern have strong effects on its 
morphology. This effect is expected to be strongest for struc- 
tures that eventually form cavities in an almost solid body. 
In fact the dark matter distribution shows much more scat- 
ter at larger radii, compared to the fairly filamentary halo 
distribution. While the scatter described so far is an intrin- 
sic feature of the point process, the very large error for the 
volume functionals at small radii is due to the Monte Carlo 
integration used for the numerical evaluation; by increas- 
ing the number of Monte Carlo points, it can be somewhat 
reduced but not fully eliminated. 

It is interesting to discuss the anti-bias of halo distribution 
observed in all of the Minkowski functionals at z = 0. Note 
that our halo finder algorithm identifies both isolated halos 
and halo s located within virial radii of other halos (satel- 
lites; see Colin et al. 1999 for details). Therefore, the small- 
scale anti-bias is not caused by the characteristic size of the 
halos and their mutual exclusion. Note, for example, that 
there is no indica tion of similar anti-bias at z = 4. Rather, 
as was shown by Kravtsov & Klypin (1999), the anti-bias 



is most likely due to dynamical friction and (to a lesser de- 
gree) to the tidal stripping in the high-density environments 



tests of the origin and evolution of halo anti-bias. Here we 
note that if the anti-bias is indeed due to the dynamical 
processes affecting galaxy evolution in groups and clusters, 
studies of higher-order small-scale galaxy clustering and its 
dependence on morphology may prove to be useful and inde- 
pendent probes of galaxy evolution in dense environments. 



5 SUMMARY AND CONCLUSIONS 

We have used a high-resolution dissipationless simulation 
of the formation of dark matter halos in a flat cosmologi- 
cal model dominated by the cosmological constant and cold 
dark matter. In the 60/i~^ Mpc simulation box we have iden- 
tified 5575, 9644, and 5869 halos of maximum circular ve- 
locity Ucirc > lllkm/s at redshifts 4, 1, and 0, respectively. 

At each of these epochs we have calculated the Minkowski 
functionals for bodies constructed as a Boolean grain model 
of 5000 randomly selected halos and 5000 randomly selected 
dark matter particles, respectively. The Minkowski function- 
als focus on global features of three-dimensional bodies, and 
provide a unifying frame for the analysis of cosmic struc- 
tures which comprises the void probability function as well 
as the genus statistics. Minkowski functionals do not replace 
but complement such traditional tools as the two-point cor- 
relation function, since they reveal clustering properties of 
higher-order. 

We found almost no evolution of the Minkowski function- 
als of halos between z = A and 2 = 1 and little evolution 
between 2 = 1 and z = 0. The rather small differences at 
various epochs are mainly due to evolution of nonlinear clus- 
tering on small scales. Already ai z — 4 the Minkowski 
functionals of halos are markedly different from a Poisson 
process. We can deduce that from the beginning, luminous 
matter quickly forms complicated structure on large scales, 
which may be interpreted as a network of sheets and fila- 
ments. 

On the other hand, the clustering of dark matter evolves 
considerably with time. At 2: = 4 the clustering properties 
of DM mostly resemble those of the Poisson process, with 
only minor deviations. At the present epoch, the clustering is 
as strong as for the halo sample, although the morphology 
of the dark matter distribution differs from the structure 
formed by the halos. 

The difference of the Minkowski functionals of halos and 
dark matter can be interpreted as a scale dependent bias 
which is qualitatively similar but quantitatively different for 
each of the Minkowski functionals. The evolution of halo bias 
in the linear and mildly nonlinear regimes results from an 
interplay between halo formation and merging rates in dif- 
ferent regions and in the field. Both rates are non-monotonic 
functions of time, but their combined effect is decrease of the 
initially high statistical bias during the course of the evolu- 
tion. In the non-linear regime, the halo bias evolution may 
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be also affec ted by the dynamical pr ocesses inside ciusters 
and groups (Kravtsov & Klypin 1999). 



We have demonstrated in this paper that the Minkowski 
functionals contain information about the correlation func- 
tions of arbitrary order. Therefore, our results imply that 
the scale-dependent bias is a property of not only the two- 
point halo correlation function, but also of the correlation 
functions of higher order. This, in turn, means that there 
are morphological differences between distributions of halos 
and matter at both small and large scales. In particular, at 
early epochs the halos are already strongly clustered and 
their distribution shows distinct features of a network. The 
dark matter distribution at the same epochs is still weakly 
clustered and is close to the Poisson distribution. At the 
present epoch the halo distribution appears less clustered 
on small scales, while it remains more filamentary on large 
scales than the distribution of matter. Whether the mea- 
sured Minkowski functionals indicate the abundance of fila- 
ments or sheets or both in the network is currently an open 
question. Although the very steep initial rise of the curves in 
Figure H would favor a filamentary distribution, sheet-like 
structures are not ruled out either. Quantitative measure- 
ments of planarity and filamentarity will be the subject of a 
follow-up article. 
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APPENDIX A: MINKOWSKI FUNCTIONALS 
EXPRESSED IN TERMS OF CORRELATION 
FUNCTIONS 

It is known that Minkowski functionals contain correlation 
functions of arbitrary order. Here we aim to establish a pre- 
cise relation between the average Minkowski functionals t;^ 
per unit volume in terms of the hierarchy of correlation func- 
tions ^n of the point process. 

In order to cut down on notation we always use a linear com- 
bination of Minko wski functionals, the so-called Minkowski 
polynomial M{t) ( [Hadwiger 1975| ; Kellerer 198^ ). We define 

d 
/j=0 



Qui^ 



V,iK) 



(Al) 



with coefBcientsH 



i^d~ 



(^d 



(A2) 



Analogously, we construct a polynomial m(t) from the den- 
sities v^ of the Minkowski functionals per unit volume. From 
now on, we drop the variable t from the polynomial, and only 
reintroduce it when extracting the individual Minkowski 
functionals via 

V,{K)^—^M{t;K). (A3) 

By taking the ensemble average and applying the ad ditivity 
relation repeatedly, we obtain the formula given by Mecke 



et al. (1994) 



E 



(-1)" 



dr„p„ (xi. 



,x„)Af (Bxi n...nBx„ 



with the integration measure defined by 

dr„ = / d'^xi . . .d'^Xn. 



(A4) 



(A5) 



If the density functions p„ were independent of position, the 
integrals c ould be performe d using the principal kinema tic 
formula by Blaschke (1936) (see also Chern fc Yien 194C| for 



the more general version employed here) . Fixing one body A 
in space, and intersecting with a moving body B, we obtain 
after integration over all possible positions of _B a simple 
factorization. 



d^a; M(AnB^) = M {A)M [B), 



VW 



* LOfi denotes the volume of the /^-dimensional unit ball. Some 
important special values are ojq = 1, t<Ji = 2, uJ2 = i", and uj^ = 



where equality holds to order f*. Unfortunately, straightfor- 
ward application of the principal kinematic formula to Equa- 
tion (Kj) is impossible in most situations since the density 
functions are in general position dependent. However, it is 
possible to use the factorization properties of the density 
functions into the correlation functions. 

Equation ( |A4| ) motivates the introduction of a hierarchy of 
"moments" P„ by weighting the n-point density function 
with the Minkowski polynomial of n intersecting balls, 

Pn = / dT„p„ (xi, . . . ,x„) M (Bxi n . . . n Bx„) . (A7) 

Writing the density functions p„ as sums over all possible 
partitions into correlation functions ^„, and using the princi- 
pal kinematic formula, it can be shown that the "cumulants" 
corresponding to the moments P„ are the S„ given by 



dT„^„ (xi, . . . , x„) M (Bxi n . . . n Bx„ 



(A8) 



Given a generating functional P[j] for the P„, a generating 
functional H[j] for the H„ is obtained through 



Setting j = — 



P[j] =exp{E[pj]). 
1, we directly arrive at the result 



(-ir 



Pn 



exp K] 



-pr= 



(A9) 



(AlO) 



Recovering the Minkowski functionals from the polynomial 
via Equation (KS), the relation in d = 3 dimensions reads 

Ho = 1 - exp {~pVo) 
vi = exp(-pl/(,) pVi 



V2 = exp (-pV^o) ( pV2 



V3 



■ exp (-pVo) pVs 



StT 2f72 



Ip'V,V2 + '^p'vt 



(All) 



which looks remarkably similar to the result for the Poisson 
process, except that the bare Minkowski functionals of balls 
V^(B) are replaced by the quantities 



V, = V^iB) + J2 



(-P)" 
-(n + 1)! 



dd id 

a;i . . . d a;. 



^„+i(o,xi,... ,x„)V^(BnBxi n...nBx,J. (A12) 



Note that Equ ation (A12) is a generalization of the result by 
White (1979) for the void probability function. This shows 
that all Minkowski functionals depend on the complete hier- 
archy of correlation functions through an asymptotic series, 
hence each order n can be expected to contribute signifi- 
cantly to the overall result. In Section this claim is backed 
by a direct comparison of numerical results. 

While it looks promising to push these anal ytical considera- 



tions forwa r d towards a hierarchical m odel (Balian & Scha- 



effer 1989a ; Balian fc Schaeffer 1989b ) or a lognormal Cox 
process (CJoles fc Jones 1991), we will stop here, since the 
series to second order is sufficient for our purposes. 
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